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Abstract 



We first review previous work on anisotropic flow at the AGS and 
SPS. Then the physics related to flow is discussed as well as the inter- 
^T) • action of flow with other non-flow measurements. From 40k RQMD 

I and 100k HIJING events predictions for anisotropic flow at RHIG are 

presented. Using the STAR detector acceptance, estimates for the res- 
■ olution obtainable with STAR are shown. We conclude that it should 

0^ , be possible to obtain good measurements for elliptic How with either 

the STAR main TPG or forward TPGs. Anisotropic flow should be 
^ ' easily one of the first results from STAR. 



3 1 Introduction 

> 

^ ■ The study of collective flow in nuclear collisions at high energies has been 

. attracting increasing attention from experimentalists |l|, |^. This is partly 

because recent progress has been made in the development of new techniques 
suitable for flow studies at high energies 0, |, H, ||, |, |^. Instead of studying 
(Px), in these new methods a Fourier expansion of the azimuthal distribution 
of particles is used in which the flrst harmonic coefficient, Vi, quantifies the 
directed flow and the second harmonic coefficient, V2, quantifies the elliptic 
fiow. In some cases Al and A2 were reported, which in modern terminology, 
are twice the square of the sub-event resolution. Using these new techniques 
anisotropic fiow has now been observed for heavy symmetric systems at both 
the AGS and SPS. 

At the AGS the E877 Collaboration pioneered the use of the Fourier ex- 
pansion method to measure vi and V2. They studied these quantities (as well 
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as V4) from a calorimeter as a function of centrality in different pseudorapid- 
ity windows [Q. Then they studied nucleons as well as pions as a function 
of pseudorapidity for different centralities 0. Using their spectrometer to 
identify particles while still obtaining the event plane from the calorimeter, 
they measured vi and V2 as a function of pt for different rapidities and cen- 
tralities |10|. They also reported (p^) as a function of rapidity [10|. In their 
latest papers they extended this study to hght nuclei |]lT], |12|. The E802 
Collaboration studied (px) for light nuclei in the target rapidity region using 
a forward hodoscope to determine the event plane . 

At the SPS NA49 first observed elliptic flow in a calorimeter study which 
reported A2 as a function of centrality [14|. WA98 reported Al as a function 



of centrality for protons and 7r+ in the target rapidity region |TB|, |TB[. They 
also studied (px) in the target rapidity region [|l^]. NA45 used silicon drift 
detectors to study vi and V2 as a function of pseudorapidity [0. NA49 has 
presented a differential study of vi and f 2 as a function of pt and y and 
has also started to study the centrality dependence 
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Also, the importance of flow for other measurements has just begun to be 
studied. For two particle correlations relative to the event plane the mathe- 

0] 



20, 21 



Some first results have 



matical scheme has been worked out |T! 

been given by WA98 [jl5l. Also, for non-identical particles the correlation 
relative to the event plane has been discussed [El . 



2 Physics Motivation 

Anisotropic flow, in particular elliptic flow, in spite of the relatively small 
absolute value of the effect, contains very rich physics. In general words, it 
is very sensitive to the equation of state which governs the evolution of the 
system created in the nuclear collision. Being such, anisotropic flow provides 
important information on the state of matter under the extreme conditions 
of the nuclear collision. The anticipated phase transition to QGP should 
have a dramatic effect on elliptic flow due to the softening of the equation of 
state. 

First it was pointed out in the pioneering work of 011itrault[|l, who sug- 
gested elliptic anisotropy as a possible signature of transverse collective flow. 
Within the hydro-dynamical model Ollitrault analyzed the role of different 
equations of state and phase transitions on the flnal anisotropy. Hung and 



Shuryak |25| suggested scanning with beam energy in order to look for the 
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QCD phase transition. Using their idea of the softest point in the equation 
of state combined with hydro-dynamical calculations, Rischke [|26| predicted 
a dramatic drop in the elliptic flow signal at the corresponding beam energies 



(in the original calculations this was at AGS energies). Sorge has shown|27 



that the elliptic flow is very sensitive to the pressure at maximum compres- 
sion, which is the most interesting time in the system evolution. Recent stud- 
ies |p8| within the parton cascade model yield similar conclusions providing 



also the relation between the strength of the elliptic flow and parton-parton 
cross sections. Recently, Sorge also tried p9| to combine the early system 
evolution in accordance to a QGP equation of state with a later hadron cas- 
cade. He looked at the centrality dependence of the elliptic flow in order 
to detect QGP production. Summarizing this part, we would conclude that 
the effect of QGP should be seen in the anisotropic flow dependence on the 
energy of the colliding nuclei, or in the dependence on the centrality of the 
collision. If the situation would be such that a QGP is produced only in a 
small fraction of the collisions than fluctuations in flow would be one of the 
best observables for this effect. 

The formation of DCC in nuclear collisions could also result in an event 



anisotropy. It could be due to the anisotropic shape of the DCC domains ||30 
or just to local fluctuations in the charged multiplicity, which should result 
in "orthogonal" flow in charged and neutral sectors . 

The very magnitude of anisotropic flow is sensitive to the degree of equi- 
libration in the system. Note that at present there is no calculation based on 
the hydro-dynamical picture which accounts for the experimentally observed 
values of the effect. This could have its origin in the obvious difficulties of hy- 
drodynamic model calculations, but it could also indicate a non-applicability 
of the picture to nuclear collisions. The cascade models such as RQMD 
describe the data much better. From this point of view the analysis of el- 
liptic flow in the collision-less and hydrodynamic limits performed in ||22 
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very interesting. The HBT interferometry performed relative to the event 
plane [1^, [2^, |21|, |2^, |23| becomes also extremely important at this point. 
Does the system really expand in the reaction plane as prescribed by hy- 
drodynamics? Simultaneous measurements of the anisotropic flow and the 
two-particle, identical as well as non identical [^], correlations in principle 
should answer this question. 

We must also mention the importance of anisotropic flow measurements 
to the vast variety of other measurements, which from first look have nothing 
to do with anisotropic flow. Let us consider high pt particle production. It 
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could be that the production mechanism (hard parton scattering) is very 
insensitive to the in-plane expansion, but that the rescattering of high pt 
partons is different in the different directions of particle emission due to the 
anisotropic geometry of the collision zone. This would lead to anisotropy in 
high Pt particle production and gives another opportunity to study how it 



develops |2|, gg. 

Another example is HBT measurements averaged over all orientations 
of particle emission. One would think that this does not require reaction 
plane measurements, but this is not really true. The mixed pair distribution 
usually used in the correlation function calculation can strongly depend on 
the relative orientation of the reaction plane of the events used to create 
the mixed pair. Therefore one should have this information even in the case 
where the dependence of the HBT parameters on the reaction plane is not 
studied. 



3 Technical Requirements 

The study of azimuthal anisotropy of unidentified charged hadrons needs the 
momenta of the particles but does not have any unusual requirements for 
calibrations, momentum resolution, acceptance, efficiency, two-track resolu- 
tion, or two-track efficiency. However, for future analyses it would be good 
to have particle identification. 



4 Directed and Elliptic Flow at RHIC 

The anisotropy in the azimuthal distribution of particles is often charac- 
terized by Vi, V2 and called directed and elliptic flow respectively. This 
anisotropy, especially V2, plays an important role in high energy nuclear col- 



lisions and is expected to be even more important at RHIC energies . The 



azimuthal distribution of particles is described by a Fourier expansion |Q 

where is the true reaction plane angle. The reaction plane is defined by 
the beam direction and the impact parameter vector b. In a given rapidity 
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(y) and pt interval the coefficients are determined by 



w„= (cos[n(0-^,)]). (2) 

Similarly this Fourier expansion can be done in coordinate space, where for 
a given rapidity and pt interval the coefficients are determined by 

y 

Tn = (cos[n(arctan(— ) — '^r)]) (3) 

where x^y are the particle space coordinates at freeze-out. Of course, these 
equations only apply to simulations where one knows 

Comparing the anisotropy coefficients in momentum space with the 
anisotropy coefficients in coordinate space (r„) as a function of pt helps us 
to understand the space-time evolution of nucleus-nucleus collisions |[T^, |3^ . 
To study this space-time evolution at RHIC, 40 000 Au+Au collisions at 

= 200 AGeV have been analyzed using the RQMD v2.4 model. 

Figs. [l|a-d show the ffist harmonic both in momentum and coordinate 
space for nucleons and pions. For nucleons at mid-rapidity note the similar- 
ity in shape of Vi versus y and ri versus y. Here (Fig. la) both the slopes of 
f 1 versus y and ri versus y show a reversal of sign. This finds an explanation 
in a picture with strong (positive) space-momentum correlations, taking into 
account the correlation between nucleon stopping and the original position 
of the nucleons in the transverse plane |^ . For pions, the rapidity depen- 



dence of Vi is predominantly governed by rescattering on comoving nucleons. 
Figs. |I|e-h show V2 for nucleons and pions. For both nucleons and pions V2 
is positive and is larger for particles with pt > 1.5 GeV. Particles acquire a 
large pt when they are produced by a hard collision (which should not pro- 
duce an event anisotropy) or when they have a large number of soft collisions 
(rescattering). The latter would explain the increase in V2 and it explains 
why T2 goes from negative for nucleons integrated over all pt to positive for 
nucleons with large pt- 

Collective flow and the coefficients vi and V2 are usually associated with 
soft processes. However, the coefficients describe the event anisotropy and 
are not limited to only soft physics. At RHIC energies hard processes become 
important. They happen early in the reaction and thus can be used to probe 
the early stage of the evolution of a dense system. During this time a quark- 
gluon plasma (QGP) could exist. Associated with hard processes are jets. 
However, when the transverse energy of the jets becomes smaller it becomes 
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Figure 1: Anisotropy coefficients 
for nucleons and charged pious in 
RQMD for collisions in the impact 
parameter range of 5 <h < 10 fm. 
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Figure 2: Anisotropy coefficients 
for nucleons and charged pions in 
HIJING for collisions in the im- 
pact parameter range of 5 < b < 
10 fm. 



increasingly difficult to resolve them from the "soft" particles. These jets 
with Et < 5 GeV are usually refered to as mini-jets. At RHIC energies it 
has been estimated that 50% of the transverse energy is produced by mini- 



jets [36 



Medium induced radiative energy loss of high pt partons (jet quenching) 
could be very different in a hadronic medium and a partonic medium. Re- 
cently it was shown that this energy loss per unit distance, dE/dx, grows 
linearly with the total length of the medium [^. For non central collisions 
the hot and dense overlap region has an almond shape. This implies different 
path lengths and therefore different energy loss for particles moving in the 
in-plane versus the out-plane direction. To study this anisotropy with respect 
to the reaction plane [52], 100 000 Au+Au collisions at y/s = 200 AGeV have 
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been generated using HIJING vl.35. 

Figs. ||a-d show vi and V2 for nucleons and charged pions. The coefficient 
f 1 shows a small negative slope around mid-rapidity for both nucleons and 
pions and this becomes more pronounced for particles with pt > 1.5 GeV. 
The coefficient f 2 is slightly negative over the whole rapidity range for both 
charged pions and nucleons. For particles with pt > 1.5 GeV, V2 becomes 
more negative especially at forward and backward rapidity. Figs. §e-f show 
that without jet quenching the anisotropy coefficients become zero. This 
indicates that interactions among particles, either quenching or rescattering, 
are important for producing the anisotropy. 



5 Event Plane Resolutions 

Within event generators the true reaction plane angle is known. This is 
not the case experimentally and the reaction plane has to be estimated from 
the data. This is done using the anisotropy in the azimuthal distribution 
of particles itself. The estimated reaction plane angle for the n*'* harmonic 
is called "^n- The magnitude of the anisotropy and the finite number of 
particles available to determine this event plane leads to a finite resolution. 
Therefore, the measured f coefficients with respect to the event plane have 
to be corrected for this event plane resolution 



" (cos[n(^„ - ^r)]) ■ ' ' 

However, eq. ^ uses the true reaction plane which is not known experimen- 
tally. Following Ref. 0, if one constructs the event plane from two random 
subevents one can relate the resolution of the subevents to the full event 
plane resolution, 

(cos[n(vl/„ - ^r)]) =Cx v/(cosKvI/^-vl/^)]), (5) 

where C is a correction for the difference in subevent multiplicity com- 
pared to the full event and \E'^, are the angles of the event planes deter- 
mined in the subevents. 

To calculate how well the event plane can be determined in STAR, we 
considered the TPC (-1.5 < ?/ < 1.5) and the FTPCs (2.5 < |y| < 4.). For this 
the RQMD v2.4 model predictions for Au+Au at y/s = 200 AGeV have been 
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Figure 3: RQMD v2.4 prediction for V2 using 
TT"*" + 7r~ within -1.5 < y < 1.5. The multiplicity 
and event plane resolution are also shown. 



used. In Fig. |a, f 2 for charged pions integrated over the TPC rapidity region 
is shown versus the impact parameter b. Fig. shows the corresponding 
multiphcity as a function of b. These quantities lead to a resolution for t>2, 
calculated using the true reaction plane, as shown in Fig. The resolution 
for V2 which can be obtained in the STAR TPC using subevents is shown 
in Fig. I^d. For V2 charged pions and protons both contribute positively and 
therefore do not need to be identified. However, the multiplicity of protons 
at mid-rapidity is small compared to that of pions and, therefore, including 
protons does not significantly change the resolution. 

In Fig. I^a, f 2 integrated over the FTPC rapidity region is shown versus the 
impact parameter b. For the FTPCs the 7r+, vr" and protons are combined. It 
was shown in Fig. that f 2 is relatively flat as a function of rapidity and its 
magnitude is therefore comparable in the FTPC and TPC regions. Fig. ^ 
shows the corresponding multiplicity as a function of b for the combined 
FTPCs. These quantities lead to a resolution for t>2, calculated using the 
true reaction plane, as shown in Fig. ^c. The resolution for f 2 which can be 
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Figure 4: RQMD v2.4 prediction 
for elliptic flow using 7r+,7r^ and 
protons within 2.5 < \y\ < 4-0. 
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Figure 5: RQMD v2.4 prediction 
for directed flow using 7r+,7r^ and 
protons within 2.5 < \y\ < 4-0. 



obtained in the STAR FTPCs using subevents is shown in Fig. If only 
one FTPC would be used this resolution would be approximately \/2 smaller. 

Using V2 the event plane can be determined, however the sign of V2 is 
not determined relative to b. This sign could be determined from V2 relative 
to \E'i. Fig. |T]c shows that around mid rapidity Vi is maximally 0.5% which 
makes extremely hard to measure. From Fig. |l|a and |l|c it is clear that the 
best region to measure vi is at forward rapidity. Fig. ^a shows vi integrated 
over the FTPC rapidity region, versus b. As for V2, the 7r+,7r^ and protons 
are combined. This decreases the magnitude of Vi because their signs are 
opposite but the FTPCs are not able to separate these particles. At large 
b the magnitude of Vi becomes ~ 1% and, although this is already hard to 
measure, also the multiplicity decreases rapidly at large b. This leads to 
negligible resolution for vi at all values of b, which is shown in Fig. |^c. 



6 Conclusion 

We have investigated the feasibility of reconstructing the event plane. Both 
Fig. ^ and Fig. ^ show that it is possible to determine the second harmonic 
event plane and calculate f 2 within STAR, assuming the RQMD predictions 
(multiplicity distribution, magnitude of V2) are correct. For V2 both the TPC 
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or the FTPCs can be used. This would initially provide a cross check and 
later combining both detectors would increase the resolution. For this study 
we only need the momenta of the charged hadrons and thus anisotropic flow 
could be one of the first results from STAR. For future analyses it would be 
good to have particle identification. Because it is important to study the 
dependence of V2 as a function of h we would like to have 10 centrality 
bins, which would be possible with 1 000 000 minimum bias events. 
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